Understanding the dynamics of fluid–structure interaction with an Air Deflected Microfluidic Chip (ADMC)

A deformable microfluidic system and a fluidic dynamic model have been successfully coupled to understand the dynamic fluid–structure interaction in transient flow, designed to understand the dentine hypersensitivity caused by hydrodynamic theory. The Polydimethylsiloxane thin sidewalls of the microfluidic chip are deformed with air pressure ranging from 50 to 500 mbar to move the liquid meniscus in the central liquid channel. The experiments show that the meniscus sharply increased in the first 10th of second and the increase is nonlinearly proportional to the applied pressure. A theoretical model is developed based on the unsteady Bernoulli equation and can well predict the ending point of the liquid displacement as well as the dynamics process, regardless of the wall thickness. Moreover, an overshooting and oscillation phenomenon is observed by reducing the head loss coefficient by a few orders which could be the key to explain the dentine hypersensitivity caused by the liquid movement in the dentine tubules.


Scientific Reports
| (2022) 12:20399 | https://doi.org/10.1038/s41598-022-24112-w www.nature.com/scientificreports/ In this work, we construct an Air Deflected Microfluidic Chip (ADMC) with PDMS and study the fluid displacement by inputting the air pressure on the deformable sidewalls. The meniscus height change is characterized by a commercial available goniometer and a custom designed LabVIEW program is used to process the collected images. The present study is designed to fill the knowledge gap in fluid-structure interaction in transient flows by performing microfluidics experiment and develop a theoretical model that can match and explain the fluid dynamics observation, thus establishing a key step for the understanding of dentine hypersensitivity in microscale.

Experiments
Negative mold fabrication. The negative mold was fabricated using an ELEGOO Saturn mono-stereolithography (MSLA) 3-D printer with Accura 25 resin (Fig. 1a-i). Due to the high hydrophobicity, the Accura 25 resin allows for demolding of cured PDMS without salinization. High aspect features such as large length versus small width/height and large height versus void geometry were achieved at a resolution of 50 µm. The complete mold was constructed out of 3-D printed sidewalls with grooves to hold 4 typical glass microscope slides. The wall structure can be freely placed on the negative mold with a loose tolerance. Use of a sealing membrane such as aluminum foil secured by tape avoids PDMS leaking during mold casting and curing. The mold and wall structure are reusable provided cleaning between each casting. In this work, we tested the samples with an aspect ratio (height/thickness) ranging from 1.5:1, 2:1, and 3:1.
Microfluidic chip fabrication. PDMS base to crosslinker with a ratio of 10:1 was used based on the instruction from the manufacturer. We casted 10 mL of PDMS into the negative mold and a standard 25 psi vacuum desiccator was used to remove air bubbles from the uncured PDMS, paying special attention to removing air bubbles in void gaps ( Fig. 1a-ii). The casted PDMS and mold were placed in a free convection oven at 45 °C for at least 4 h. Low temperature is required to avoid entering glass transition and deformation of the negative mold. After demolding, a biopsy punch (diameter: 1.5 mm) was used to form inlets and outlets in the desired channels ( Fig. 1a-iii). For chip completion, we used a spin coated PDMS layer of approximately 1 mm in thickness as a bottom layer of the chip. To achieve this 1 mm layer, we spin coated four layers of PDMS at 350 RPM for 1 min per layer, curing in between each layer to achieve a layer of PDMS with an approximate thickness of 1 mm. Then, a very thin film of PDMS was spin coated on the cured 1 mm layer at 1000 RPM for 1 min as an adhesive layer between the spin coated PDMS and the casted PDMS chip. Finally, a glass capillary and an inlet plug were inserted into the ADMC chip ( Fig. 1a-iv). A photograph of the dye filled ADMC sample is shown in Fig. 1b Measurements of applied pressure and fluid displacement. The experimental setup is shown in Fig. 1c. Reservoir pressure was supplied to the system using lab air connected through a Drierite gas purifier (desiccator) inline to an Elve Flow flow controller (OB1 MK3+ to keep the air supply dry and clean. The flow controller operates between 0 and 2000 mbar and delivers reservoir pressure with a resolution of 100 µbar to the gas channels on the ADMC through rigid pneumatic tubing. The flow controller was outfitted with pressure transducers to measure applied pressure. The liquid meniscus height in the glass capillary tube was measured using a 5X-120X microscope lens attached to a 100 FPS goniometer camera. The acquisition sampling time was set to 60-75 Hz. A LabVIEW program was designed to synchronize the frame collection with the pressure measurements. Before the experiments, the flow controller and the goniometer were calibrated. Detecting fluid displacement by image processing. An interactive routine was created to streamline and control processing of the image frames and pressure data. MATLAB was used to track the position of the liquid meniscus by means of image processing (the MATLAB code available in GITHUB). Each image was converted into an edge image using a Canny edge detection function, where the contrasting edges were converted to white lines on a black background. A sensitivity threshold was set between 0.2 and 0.6 and tuned depending on light intensity. To convert the image to engineering units, we performed a centroid calculation to determine the displacement of the meniscus of that frame relative to the first. A pixel conversion factor was calculated for each dataset based on the known thickness of the glass capillary in observation. Each glass capillary has an outside diameter of 2 mm ± 0.1 mm. Compiling each of the frame data points together yielded the fluid dynamic response of our chip.
Theoretical model. A theoretical model was developed to understand the dynamics of fluid-structure interaction and to verify the experimental results. The model ( Fig. 1d-I,d-ii) generally includes two sectionsthe pressurized channel and the vertical tube at the end of the channel. Applying the unsteady Bernoulli equation to the device, where P is the pressure, ρ is the fluid density, µ is the velocity, t is time, z is the height of the fluid (z I = 0 and z T is the height of the free surface in the vertical tube), and h L is the head loss in the flow, which generally takes account of the major and minor losses in the device. The subscription T and I denote the variables in the vertical tube and the horizontal pressurized channel. The total fluid in the device is conserved, therefore,  where , and assume K = K I = K T . This system of equations was solved using the ode15i function in MATLAB. The initial velocity was assumed zero and no deformation was assumed at t = 0. Figure 2a shows the raw images labelled with 1, 2, 3, and 4 correspond images captured at 0, 0.10, 2.34, and 23.77 s, respectively. Height values that were calculated are shown alongside the time values. The frames from Fig. 2a were processed further using a Canny edge detection algorithm. In addition to the canny edge detection, the images were cropped to omit extra and unneeded data, as shown by the white rectangle in the frames. Figure 2c shows the cropped areas from Fig. 2b. The top and bottom edges of the liquid meniscus were averaged together by centroid calculation to output the liquid surface height relative to the first frame. The average meniscus height is shown in magenta in Fig. 2c. Each image in Fig. 2a-c correspond to the data points marked in Fig. 2d. Figure 2d shows a sample data set from a 500 μm wall thickness device with 200 mbar of pressure applied. The measured meniscus height versus time under various air pressure is shown in Fig. 3. For 500 µm thick sample ( Fig. 3a-I,a-ii), the meniscus shows a sharp increase within the first 10th of s, regardless of the input air pressure. In the second phase, the meniscus gradually rises over the course of ~ 0.5 s. Then, the meniscus reaches to a steady state without showing much change from 1 to 10 s. With 500 mbar air pressure, the meniscus height reaches to ~ 12,000 µm in contrast to 1000 µm for 50 mbar. As shown in Fig. 3b,c, the 750 µm and 1000 µm samples follow a similar trend with the 500 µm sample. However, the maximum meniscus height only reaches 6,500 and 1700 for 750 µm and 1000 µm samples, respectively.

Results and discussions
The experimental results and simulation comparison for 500, 750, and 1000 µm samples are shown in Fig. 4a-c, respectively. Our theoretical model well captures the dynamics of the second and third phases of the dynamic process, regardless of the wall thickness. Specifically, the heights at the end of Phase 1 were used as the starting points of the modeling and the other used parameters in the theoretical model to compare with the experimental results are listed in Table 1. The theoretical model does not only compare well with the ending points of the experimental results but also with the dynamics. The change of the rising rate and the magnitudes are well reproduced. This indicates that the theoretical model can, to an extent, reflect the complicated dynamics in the experiments in addition to the steady state results that are dictated by the mass conservation before and after the deformation.
(2) www.nature.com/scientificreports/ Much effort was also made to capture the dynamics of Phase 1, but it turns out that the meniscus height rose too quickly to compare with the theoretical model. This indicates that the process in Phase 1 was governed by a different mechanism beyond the framework of the hydrodynamics described in the theoretical model. We suspect this process was formed due to the pre-tension during the fabrication process-when the pressure was applied, the pre-deformation played a significant role to accelerate the meniscus' rising. A further investigation with different fabrication stress is needed to resolve this issue.
The last process of the experiment still shows a slow increasing trend, which is not captured by the theoretical model. This difference might be due to the response time of the structure deformation, while in the theoretical model the deformation was assumed to be instantaneous. In addition, the streamwise non-uniformity might be attributed to the difference.
An observation can be made in the parameters used in the theoretical model. The tube law's α and P c increases with the channel wall thickness. This is expected as the thicker wall leads to less flexible cross-sectional changes. The head loss coefficient's change is not expected as the wall thickness should not change the roughness of the channel surface. This result indicates that the minor loss is dominating the head loss.
As mentioned earlier, the present theoretical model assumes that (1) the channel deforms uniformly; (2) surface tension is neglectable at the meniscus in the tube; (3) the minor loss at the connection between the www.nature.com/scientificreports/ channel and the tube could be ignored; (4) the tube was assumed vertical; (5) the deformation of the channel was assumed to be instantaneous; and (6) the total fluid mass was assumed conserved. These assumptions may contribute to the overall uncertainty of the theoretical model's performance and future studies should be targeted to examine these over-simplifications. Interestingly, if the hydraulic loss factor is reduced by a few orders to allow stronger inertia in the flow dynamics, we observe an overshooting and oscillation process (Fig. 5). This is consistent with the 3-D computational fluidic dynamic studies by other groups 12,24 . The oscillation could be the key to the tooth induced by thermal triggers and our theoretical model well predicts its presence under the current framework. A further comparison with the existing numerical model is required to investigate the complicated mechanisms revealed in this study and this could lead to resolve the mystery of the particle-contained tooth paste-particle-containing tooth paste was observed effective in reducing dentine-hypersensitivity 25 . However, this mechanism has not been fully understood. The particles contained in the toothpaste could block the open tubules in damaged teeth to reduce the oscillation magnitude of the flow to reduce the pain. This hypothesis has not been tested in simulation or experiments. Our study provides a theoretical framework and experiment setup to investigate the impact of the boundary condition of the tubule flow on the flow dynamics and thus shed lights on this issue. A further study is planned to move forward in this direction. Another knowledge gap revealed in this study is the scale of the model. As the head loss coefficient could determine the presence of the oscillation of the flow and the scale of  26 , a detailed study to determine the parameters in the lab scale and the microtubule scale is required to fully address the model scaling problem. To facilitate in vivo study, microfabrication 27 and high resolution microstereolithography 28 can be used in the future to make the ADMC with a similar size of the actual dentine tubules.   Figure 5. Observed liquid meniscus oscillation at the first 0.008 s in our theoretical model.